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25 Abstract 
O 

We compute the scattering amplitude for pion scattering in Wilson chiral perturbation theory 
for two degenerate quark flavors. We consider two different regimes where the quark mass m is 
of order (i) aAq CD and (ii) a 2 AQ CD . Analytic expressions for the scattering lengths in all three 
isospin channels are given. As a result of the 0(a 2 ) terms the 1 = and 1 = 2 scattering lengths 
do not vanish in the chiral limit. Moreover, additional chiral logarithms proportional to a 2 lnM 2 
r— i are present in the one-loop results for regime (ii). These contributions significantly modify the 

familiar results from continuum chiral perturbation theory. 

^ PACS numbers: ll.15.Ha, 12.39.Fe, 12.38.Gc 
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I. INTRODUCTION 



The scattering of low-energy pions is an important process in QCD. Since the pions are 
the pseudo Goldstone bosons associated with the spontaneous breaking of chiral symmetry, 
one can use chiral perturbation theory (ChPT) to compute the scattering amplitude and the 
associated scattering lengths to high precision. These theoretical predictions are in excellent 
agreement with experiment [H |2] . 

Since a few years ago Lattice QCD has started to impact ChPT and the study of pion scat- 
tering. Increase in computer power and algorithmic advances nowadays allow unquenched 
simulations with light pion masses so that contact with ChPT can be made. Monitoring the 
quark mass dependence of the pion mass and the pion decay constant various lattice groups 
have obtained estimates for the low energy constants 1% and U that enter the one-loop ChPT 
results of these two observables [31 SI E] - These two constants play an important role for 
pion scattering, in particular Z 3 , since it dominates the uncertainties in the isospin zero and 
isospin two scattering lengths [6] (see also [2]). 

The scattering lengths can also be computed directly on the lattice, and two unquenched 
calculations for the 1 = 2 scattering length a 2 , have been performed so far. The CP-PACS 
collaboration simulated two dynamical Wilson quarks at three lattice spacings between 0.1 
and 0.2 fm [7]. The pion masses, however, were rather heavy, spanning the range 0.5 to 1.1 
GeV. Thus, contact to ChPT is not expected. The NPLQCD collaboration [91 [10| performed 
a mixed action simulation with domain- wall valence quarks and staggered sea quarks. 1 The 
four pion masses were in the range 290 to 490 MeV and the lattice spacing a ~ 0.125 
fm. 2 Since data at one lattice spacing only was generated and due to the use of different 
actions in the sea and valence sector mixed action ChPT [TH [121 H3] was used to extrapolate 
simultaneously to the continuum and chiral limit. The obtained result for a 2 , agrees very 
well with the experimentally measured one. 

In this paper we study pion scattering in Wilson ChPT (WChPT) [Ml US], the low-energy 
effective theory for Lattice QCD with Wilson quarks. We compute the scattering amplitude 
to one loop including the leading corrections due to the non-zero lattice spacing. From 
the scattering amplitude we derive expressions for the scattering lengths in all three isospin 
channels which can be used to analyze lattice data before the continuum limit has been 
taken. We work in two different quark mass regimes: (i) m ~ aAq CD and (ii) m ~ a 2 AQ CD . 
These are most likely applicable to todays and future simulations with light Wilson quarks. 
The 1 = 2 scattering length in regime (i) has already been calculated before in Ref. [T6] . 
and we agree with this result. 

The lattice spacing is, just as the quark mass, a source of explicit chiral symmetry break- 
ing. Therefore, the lattice artifacts modify the results from continuum ChPT. For example, 



1 The unquenched configurations were generated by the MILC collaboration |8] and are publicly available. 

2 The pion masses refer to the staggered Goldstone pion. The masses of the other pions are significantly 
heavier. The taste singlet pion, which is the heaviest, is about 450 MeV heavier than the Goldstone pion 
at this lattice spacing pj]. 
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the scattering lengths Oq for / = 0, 2 do not vanish in the chiral limit but assume non-zero 
values. This is not unexpected [UJ [TBI HH] - However, the value in the chiral limit is of order 
a 2 and not of order a, as one might naively think based on the (broken) symmetries of the 
Wilson fermion action. 

In addition, at sufficiently small quark masses there appear additional chiral logarithms 
proportional to a 2 In M 2 in the chiral expansion. These additional contributions may obscure 
the continuum chiral logarithms. Consequently, fits of the continuum ChPT expressions 
to the lattice data may easily result in erroneous determinations of the Gasser-Leutwyler 
coefficients associated with pion scattering. 

II. PION SCATTERING AT TREE LEVEL 
A. Setup 

The chiral effective Lagrangian of WChPT is expanded in powers of (small) pion momenta 
p 2 , quark masses m and the lattice spacing a. Based on the symmetries of the underlying 
Symanzik action [201 [21] the chiral Lagrangian including all terms of 0(p A ,p 2 m, m 2 ,p 2 a, ma) 
is given in Ref. [15]. The 0(a 2 ) contributions are constructed in Ref. [22] and, independently, 
in Ref. [23] for the two-flavor case. 

In the following we will restrict ourselves to Nf = 2 with degenerate quark mass m. In 
this case the chiral Lagrangian (in Euclidean space-time) including the p 2 ,m and a 2 terms 
is found to be [221 [23] 

C 2 = ^£^£t) _ t M z + S t) _ 1 (2 ^ + W' 8 )a 2 (£ + St) 2 , (1) 

where the angled brackets denote traces over the flavor indices. The field S containing the 
pion fields is defined as usual, 

S(x) = exp (j*{z)j , n(x) = 7r a (z)y, (2) 

with the standard Pauli matrices a a . The quark mass and the lattice spacing enter through 
the combinations 

m = 2Bm, a = 2W a . (3) 

The coefficients B and / are the familiar leading order (LO) low-energy coefficients (LECs) 
from continuum ChPT [2H [23] > while Wq, Wq, W% are LECs associated with the non-zero 
lattice spacing artifacts [T31 [22] . 

Note that the mass parameter m in eq. (|3| is the so-called shifted mass [13] ■ Besides the 
dominant additive mass renormalization proportional to 1/a it also contains the correction 
of 0(a). Consequently, there is no term linear in the lattice spacing present in the chiral 
Lagrangian in eq. Q. 
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Note also that we keep the 0(a 2 ) correction in £ 2 and promote it to a LO term in the chiral 
expansion. This is justified (and necessary) for small enough quark masses such that the 
O(m) and the 0(a 2 ) terms are of the same order of magnitude, i.e. for m ~ a 2 Aq CD . We call 
this scenario the large cut-off effects (LCE) regime, in contrast to the so-called generic small 
quark mass (GSM) regime [26j EZ]; which assumes m ~ aAQ CD . Nevertheless, even though 
we almost exclusively work in the LCE regime we will be able to obtain the corresponding 
results for the GSM regime as well by appropriately expanding our final results (see section 

nig . 

It will be useful to introduce [2B] 

c 2 = -32(2^ + 1^)-^ (4) 

for the combination of LECs in front of the 0(a 2 ) term in the chiral Lagrangian. The sign of 
c 2 determines the phase diagram of the theory [H]. 3 For c 2 > there exists a second-order 
phase transition separating a phase where parity and flavor are spontaneously broken [29J. 
The charged pions are massless in this phase due to the spontaneous breaking of the flavor 
symmetry. Outside this phase the pion mass is given by 

M 2 = 2Bm-2c 2 a 2 , (5) 

and it vanishes at m = C2a 2 /B. For even smaller values of m the charged pions remain 
massless, while the neutral pion becomes massive again [T4] . 

Negative values of c 2 , on the other hand, imply a first order phase transition with a 
minimal non-vanishing pion mass. All three pions are massive for all quark masses, and the 
pion mass assumes its minimal value at m = 0, resulting in 

M 2 min = 2|c 2 |a 2 . (6) 

The magnitude and the sign of c 2 depend on the details of the underlying lattice theory, i.e. 
what particular lattice action has been chosen. However, it is not a simple task to measure c 2 
numerically. Adding a twisted mass term [30J to the theory the pion mass splitting between 
the neutral and the charged pion is equal to 2c 2 a 2 at LO in the chiral expansion (28J I3T] . 
This has been exploited by the ETM collaboration [32| [33] and an estimate for c 2 has been 
obtained. Within errors c 2 is negative if the standard Wilson fermion action and the tree- 
level Symanzik improved gauge action is used. However, the statistical uncertainties for c 2 
were rather large due to the presence of disconnected diagrams in the calculation of the 
neutral pion mass. 

B. Scattering amplitude and scattering lengths 

We are interested in the two-pion scattering process 

7i a (p) + 7r /3 (/e) — -> ttV) + n s (k'). (7) 



3 Note that our definition for C2 differs by a factor f 2 a 2 from the one in Ref. [13] . 
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This process is described by the scattering amplitude A. It is convenient to use the three 
Mandelstam variables as arguments for it, A = A(s,t,u). The scattering amplitude is 
straightforwardly calculated as the residue of the four pion pole in the four-point function. 
Starting from eq. Q we obtain the tree-level result 4 

A(s,t,u) = jj(s - Ml - 2c 2 a 2 ) . (8) 

Setting the lattice spacing to zero we recover, as expected, the familiar result of continuum 
ChPT [21]. 

Having calculated the scattering amplitude we perform the standard partial wave expan- 
sion and obtain the scattering lengths [21]. Most interesting are the scattering lengths (Zq 
for definite isospin 1 = and 1 = 2: 

32^( M »-7 2 ^ 2 )' < 9 > 

< = -T^( M » + 2C2 ° 2 )- < 10 > 

Again, for a = we recover the tree-level continuum results, first obtained by Weinberg 
|34j . For a non-zero lattice spacing, however, the continuum results are modified in such a 
way that the scattering lengths no longer vanish in the chiral limit. Instead, they assume 
non-zero values of 0(a 2 ). In other words, the ratio a^/M 2 is no longer a constant but has 
the functional form 

_ ^oo , A i n ,N 

n it 

with A[ Q being a constant and Aq being of order a 2 . Hence, the ratio a^/M 2 diverges in 
the chiral limit. This divergence has been anticipated first by Kawamoto and Smit |17j . 
However, note that the coefficient Aq is of order a 2 rather than of order a. This holds even 
for unimproved Wilson fermions, in contrast to earlier expectations [T8| [19]. 

On the other hand, the divergence in the chiral limit will only be present if c 2 > 0, 
because only in this case can the pion indeed become massless. For the opposite sign the 
pion mass cannot be smaller than the minimal value quoted in eq. ([6]), resulting in the 
following minimal values for the scattering lengths: 

12 

I I 2 2 / \ 

a 0,min = f 2 ' a 0,min = 0- (12) 

Figure [T] sketches the pion mass dependence of the scattering lengths for the two possible 
signs of C2- It seems feasible that measurements of the scattering lengths will allow more 
precise determinations of C2- Extrapolating the data for Oq to the chiral limit one may 



Recall that the Lagrangian in Eq. ([lj is given in Euclidean space-time, so one has to Wick-rotate back to 
Minkowski space in order to get the physical scattering amplitude. 
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FIG. 1: Sketch of the scattering lengths as a function of the pion mass at non-zero lattice spacing. 
The left panel shows I67r/ 2 aQ as a function of the tree level pion mass Mq. For C2 < the pion 
mass cannot be smaller than MQ min in eq. (jjjj). Nevertheless, extrapolating to the massless point 
the scattering length assumes the value — 2c2a 2 , as indicated by the dashed line. For C2 > the 
pion mass can be taken zero. At this mass the scattering length also assumes the value — 2c2a 2 , 
now with the opposite sign. The right panel shows the analogous sketch for the 1 = scattering 
length. 



directly read off C2 as the value at vanishing pion mass, even for the C2 < case. A 



practical advantage is that the calculation of ai does not involve disconnected diagrams 



which introduce large statistical uncertainties. 

The 1 = 1 channel is somewhat special. The scattering amplitude for this isospin channel 
is given by A(t, s,u) — A(u,t,s). The c 2 contribution drops out in this difference and the 
scattering length a\ is is given by 

1 Ml (13) 



2 



1 24tt/ 

This is exactly the tree level result of continuum ChPT [21] and it suggests that the scaling 
violations in a\ are very small. 5 



III. PION SCATTERING AT ONE LOOP 
A. Power counting and higher order terms 

Having promoted the 0(a 2 ) correction to the LO Lagrangian, it will contribute to the one- 
loop results of various quantities in a non-trivial way. Expanding the 0(a 2 ) term in terms 



5 Note that our definition of the scattering length a\ differs by a factor of Mq from the one in Ref. [23] . 
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of pion fields we obtain vertices proportional to 2c 2 a 2 . As part of one-loop diagrams these 
lead to non-analytic corrections proportional to 2c2a 2 \nM^, for example in the expression 
of the pion mass (see figure [2]). 



The presence of additional chiral logarithms of order a 2 In Mq is well known in staggered 
ChPT (SChPT) [3S1 E31 07J, and has been pointed out first in Ref. [23] for WChPT. They 
are considered to be one reason why the chiral logarithms known from continuum ChPT are 
not reproduced in the lattice data: The additional chiral logs obscure the non-analytic quark 
mass dependence due to the continuum chiral logs, and the naively expected behaviour is 
lost. 

The power counting in WChPT is slightly non-trivial if we take the 0(a 2 ) term at LO. The 
LO Lagrangian consists of the terms of 0(p 2 , m,a 2 ). In order to renormalize the divergencies 
of the loop diagrams we need higher order counterterms in the chiral Lagrangian. These are, 
besides the standard ones in £ 4 of continuum ChPT [2U [25] , the terms of order p 2 a 2 , ma 2 , a 4 . 
However, terms of order p 2 a, ma, a 3 are also present and formally of lower order. Hence these 
should also be included, even though their LECs do not get renormalized at one loop. To 
conclude, one-loop calculations in the LCE regime require the following terms in the chiral 
Lagrangian: 

LCE regime: T n 2 2 
LO : p , m, a 

p 2 a, ma, a 3 (14) 



NLO 
NNLO 



4 2 2 2 2 2 4 

p , p m, m , pa, ma , a 



The standard NLO terms of continuum ChPT appear here at NNLO, a consequence of the 

fact that m ~ a 2 Aq CD in the LCE regime. 6 In the GSM regime, however, we recover the 

standard ordering. Since m ~ aAQ CD > a 2 AQ CD in this regime, a one-loop calculation 

requires the following terms: 

GSM regime: T 2 
LO : p , m 

NLO : p 4 , p 2 m, m 2 , p 2 a, ma, a 



4 ^,2 ^„ „i (15) 



Note that all terms in (15) are also present in (14), even though reshuffled. This already 
indicates, that a result obtained in the LCE regime yields the corresponding result in the 
GSM regime if appropriately expanded. 



Most of the necessary terms in (14) have already been constructed. So far unknown are 
the contributions of order p 2 a 2 , ma 2 , a 3 , a 4 . Performing a standard spurion analysis with the 
spurion fields in Ref. [22] it is straightforward to construct these missing terms, at least the 
ones we need as counterterms for the observables we are interested in here, the pion mass 
and the scattering amplitude. Some details concerning the construction of these terms are 
collected in appendix [A] 

Finally, if the underlying lattice theory is non-perturbatively 0(a) improved according 
to the Symanzik improvement program [201 EI], the terms of order p 2 a,ma are absent in 

6 The terms of 0(p 2 a, ma, a 3 ) are not present in SChPT due to the axial U(l) symmetry. Consequently, 



the terms of the third row in (14) are the NLO terms 
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OC 2C2<2 2 



FIG. 2: Additional one-loop diagramm contributing to the self energy. The vertex is proportional 



C2a 2 , leading to a chiral logarithm c 2 a 2 In Mq. 



(14) and (15). In the following we will always keep these contributions in our calculation; 



the 0(a) improved result is simply obtained by dropping the appropriate terms linear in the 
lattice spacing. 



B. Pion mass to one loop 



The modification due to the additional chiral logarithms is illustrated best in the one- 
loop result for the pion mass, which we will also need in the next section. In terms of the 
tree- level pion mass Mq, cf. Eq. we find 



Mi 



Ml 



1 + 



M ° 2 ln 



MS 
AS 

, , 2c 2 W a 3 (2c 2 a 2 ) 2 
+ h h h 



32tt 2 / 



+ 



5 2c 2 a 2 
32tt 2 f 2 



In 





W a~ 


(f) 





p 



p 



(16) 



The coefficients A 2 , S 2 and k\, k%, k^ are (combinations of) unknown LECs. A 2 is the familiar 
scale independent LEC present in the continuum result [21], and in complete analogy we in- 
troduced the new LEC S3. Note that we have chosen the coefficients k{ to be dimensionless. 7 
In the 0(a) improved theory the k% term is absent. 



A few comments concerning (16) are in order. First, setting the lattice spacing to zero we 
recover the familiar result from continuum ChPT [24J , as expected. This result gets modified 
at non-zero lattice spacing. In particular, there exists the anticipated chiral logarithm 
proportional to c 2 a 2 In M5 2 . Note that the coefficient in front of it is ten times larger than 
the coefficient in front of the continuum chiral log proportional to MqIuM^. Hence, even 
small 0(a 2 ) contributions may dilute the continuum chiral logarithm completely. This is 



better seen if we rewrite the square bracket in (|16|) as 
1 



1 + 



32tt : 



: {M 2 + 10c 2 a 2 }ln 



Ml 
A3 



+ 



10 



r c 2 a 



2 In 



(17) 



327T 2 / 2 ' V S 3/ P 

so that the quark mass dependence comes entirely from the fiiM^/A! term. Negative values 
of c 2 can render the factor M5 2 + 10c 2 a 2 exceptionally small such that the chiral logarithm 
is effectively not active. 



7 Details about how the LECs in (16 1 are related to the ones in the chiral Lagrangian can be found in 
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This scenario is not as unlikely as one may think. We remark that the ETM collaboration 
has found a negative value for c 2 in their twisted mass simulations, and the calculation of 
the pion mass splitting provides a rough estimate for — 2c 2 a 2 = M 2 ± — M 2 . The data [33] 
for a w 0.086 fm and M^± « 300 MeV results in -2c 2 a 2 « (185 MeV) 2 . Such a value 
completely suppresses the chiral log for pion masses around 400 MeV, a value not unusual 
in lattice simulations performed these days. 

Finally, note that the pion mass (16) does not vanish in the limit Mq = because 



of the corrections proportional to k$, k&. However, these are contributions of 0(a 3 , a 4 ) to 
the additive mass renormalization (the critical quark mass) and can be absorbed by an 
appropriate finite renormalization. To be specific, define fh and Mq by 



Mo 



2Bfh 



2Bm - 2c 2 a 2 + h 



2c 2 W a 3 
' P 



+ 



, (2c 2 a 2 ) 2 Woa t 2c 2 W a 3 \ , 
h t h^rh t — I > ( 18 ) 



P 



P 



P 



such that the quark masses m and m differ by order a 2 and higher. In terms of fh the result 



(16) reads (up to the order we are working here) 



Mi 



Ml 



1 



32vr 2 p 



Ml 



2c 2 a 2 



32tt 2 p 



In 



M| 



+ h 



W a 
~P 



(19) 



The ^3, k 4 contribution no longer appear explicitly but are absorbed in the definition of the 
quark mass fh. With this parametrization the pion mass vanishes in the chiral limit for 
rh — 0. Therefore, fh is proportional to the subtracted bare lattice quark mass (mo — m cr ) 
if the critical quark mass m cr is defined by a vanishing pion mass for m = m CT . In other 



words, doing the definition (18) we have appropriately matched the chiral effective theory 



to the underlying lattice theory for this particular renormalization condition. Obviously, 
the matching differs for other definitions of the lattice quark mass, for example through the 
PCAC relation. 

In the following we are not interested in the quark mass dependence of scattering observ- 
ables, but rather in the dependence on the pion mass. For this the parametrization in terms 



of m is sufficient, since at the end we will replace m by M 2 using Eq. (16). 



C. Scattering at one loop 

At one loop there contribute the four diagrams in figure [3] to the four-point function. 
The vertices in these diagrams can be either a vertex also present in continuum ChPT or 
the vertex proportional to c 2 a 2 stemming from the 0(a 2 ) term in Eq. (Jl|. Hence, we expect 
additional chiral logarithms c 2 a 2 MQ In Mq and (c 2 a 2 ) 2 In Mq in the scattering amplitude and 
the scattering lengths. 

The one-loop result for the scattering amplitude can be written as 

A(s,t,u) = -?-(s-M Q 2 -2c 2 a 2 ) + B(s,t,u) + C(s,t,u), (20) 
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FIG. 3: One-loop diagramms contributing to the four-point function. The diagrams with one or 
two vertices stemming from the 0(a 2 ) term in the chiral Lagrangian will give rise to chiral logs 
proportional 02a 2 Mq lnMp and (c2d 2 ) 2 In Mq . 



where the functions B and C contain the one-loop corrections. The general structure of 
these functions is 

B(s,t,u) = B COQt (s,t,u) + 2c 2 a 2 B a 2(s,t,u) , (21) 
C(s,t,u) = C cont (s,t,u) + 2c 2 a 2 C a 2(s,t,u) . (22) 

The first parts B cont , C cont are the contributions from continuum ChPT and can be found 
in Refs. [2H 139]. The contributions B a 2,C a 2 come from the 0(a 2 ) correction. Since these 
expressions are cumbersome and not very illuminating we present them in appendix [B] Here 
we just mention that these functions are finite in the limit a — » 0, so in the continuum limit 
we recover the continuum result for the scattering amplitude. 

With the amplitude at hand we can compute the scattering lengths as before. As a 
short-hand notation we introduce 

k(M 2 ) = -In (^pj , i = 1,2,3 (23) 

&(M 2 ) = -ln(|j^ , i = 3,. ..,6 (24) 

for the chiral logarithms in the following expressions. The U are the standard scale invariant 
LECs of continuum ChPT [21], and the parameters & are defined in complete analogy to 
them. However, we have made the dependence on Mq explicit. This seems appropriate since 
the pion mass is varied in lattice simulations and the h,^i are indeed functions of the pion 
mass. 

With these definitions the one-loop results for the isospin zero and isospin two scattering 
lengths read 

7 M 2 ( 5 Ml f=- .0. 21 



32tt P 



(l+8^^{'^) + 2^) + |} 



5 2c 2 a 2 / W a 



i-^-^{nea(M 2 )-i}) 



10 



Ml 



M 2 ), 



(25) 



16tt/ 2 
2c 2 a 2 
16tt/ 2 



1 - 



Ml 

127T 2 / 2 

W a 



P 



-fa 



h(Ml) + 2l 2 (Ml) + 

{llUMD-7} 



2c 2 a 2 

327T 2 / 2 



{ll£ 5 (M 2 ) + 2} 



2c 9 a 2 



32vr 2 / 2 



(26) 



Here fa is another LEC associated with the 0(a 3 ) terms in the chiral Lagrangian. As 
anticipated, we find additional chiral logarithms a 2 Ml In Mq (the terms involving £4, £ 5 ) and 
a 4 lnM 2 (the £ 6 term). 

The results (25) and (26) are given as a function of the quark mass m via the tree- level 
pion mass Mq, except for the last line in Oq, where the one-loop expression M 2 enters. 8 As 
already mentioned at the end of section IIIB[ these results have to be properly matched 
to the lattice theory. Depending on the particular renormalization condition for the quark 
mass the final results differ by terms of order a. For this reason it is advantageous to express 
the scattering lengths as a function of the pion mass and not of the quark mass. Using the 
one-loop result (16) in order to replace M5 2 by M 2 we obtain the scattering lengths as a 
function of M 2 : 



J_Ml 
32vr P 



5 2c 2 a 2 

32tT f 2 

Ml 



1 + 



5 Ml 
84vr 2 P 



5 W a 



h(M 2 ) + 2/ 2 (M 2 ) - lUMl) + 



21 



5 2c 2 a 2 f 15 



84vr 2 p 



1 



ST, , U\ W ° a U 
{fa + fa\ — fa 



2c 2 a 2 



UmD + vumI 

2c 2 a 2 



y e 5 (M 2 ) + — 



p 



P 32n 2 P 



{lie 6 (M 2 )-l} 



(27) 



16tt/ 2 



2c 2 a 2 
16tt/ 2 



Ml 



fa 



12tt 2 P 
W a 



UMl) + 2UMl)-lUMl) + l 



P 



{fa + fa} 



2c 2 a 2 

327T 2 / 2 

W a 



{^UMl) + llUMl) + 2} 



P 



k 



2c 2 a 2 

: P 32tt 2 / 



2c 2 a 2 

2 f2 



{uUmD-7} 



(28) 



Here the coefficients k and ^ are as in (23), (24) but with the tree- level pion mass replaced 
by M 2 . 

Note that we have still expressed the scattering lengths as a function of /, the decay 
constant in the chiral and continuum limit, and not in terms of f^. The reason is that f n 
has not been computed yet to one loop in the LCE regime. 9 This is not an impediment to 



8 The definition of et[] involves the physical pion mass when one goes on-shell and uses s + t + u = 4M% . 

9 The calculation of the decay constant poses additional complications: The axial vector current in the 
effective theory needs to be constructed to the appropriate order and the proper renormalization condition 
has to be taken into account [10] . 
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our results, but it means that / is a free fit parameter when our results are used to fit them 



to numerical lattice data. 10 We briefly come back to this in section HI E 



Our final results for the scattering lengths differ significantly from the corresponding ones 
in continuum ChPT. The scattering lengths do not vanish in the chiral limit, but rather 
diverge as a 4 lnM 2 . However, if C2 < there is no divergence since the pion mass cannot 
become smaller than the minimal value in eq. ([6]). For the opposite sign massless pions 
are in principle possible, but for small pion masses the chiral expansion eventually breaks 
down when C2a 2 lnM 2 becomes of order unity. In that case higher order terms leading to 
powers (c2a 2 In M 2 ) n , n = 2,3, . . . become relevant too and a summation of all these terms 
is necessary. This can presumably be done along the lines in Ref. (23], where a resummed 
pion mass formula has been derived. Here, however, we assume that the pion masses are 
heavy enough such that the chiral logarithm C2a 2 lnM 2 is a reasonably small correction to 
the leading order contribution. This is most likely the relevant case for actual numerical 
simulations. 

In any case, the pion mass dependence of the scattering lengths on the lattice can be 
very different in contrast to what one may expect from continuum ChPT. Consequently, 
attempts to fit lattice data using expressions from continuum ChPT may easily fail. 



D. Results for the GSM regime 

In this section we summarize the results for the GSM regime where m ~ aAq CD > 
a 2 AQ CD . These are easily obtained from the expressions in the last section by expanding 
the logarithms according to Ii^Mq) « \n(2Bm) — 2c2a 2 /2Bm and dropping consistently all 
higher order terms. 



For example, the one-loop expression (|16|) for the pion mass reduces to 

- 2c 2 a 2 . (29) 



Mi = 2Bm 



2Bm , /25m \ , W a 



32tt 2 / 2 V A 3 / P . 



This result is easily understood. The LO Lagrangian consists only of the kinetic and the 
mass term. Hence, the tree-level pion mass is equal to 2Bm and the chiral logarithm in the 
one-loop result is just the one from continuum ChPT. The 0(ma,a 2 ) terms enter at NLO 
and give analytic corrections only. 

Similarly we obtain the results for the scattering lengths: 



5 2c 2 a 
32vr f 2 



2 



(30) 



10 Note that the replacement of / by / w does not reduce the number of unknown LECs since the latter 
contains new LECs that are not present in the expressions paramterized by / |40j . 
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2c 2 a 2 



16np 



W a\ 

l— FT- 1 

(31) 



As expected, the results reduce to the continuum one-loop result plus analytic corrections 
of order M 2 a and a 2 . 



E. Practical remarks 

An important application of WChPT is to provide formulae that can be used to fit lattice 
data and thereby allow the chiral extrapolation to the physical pion mass. As a by-product 
one also obtains estimates for the Gasser-Leutwyler coefficients involved in these formulae. 
The results for the scattering lengths of the previous sections serve exactly this purpose, 
but the way we have written them is not practical. In order to compare and cross check our 
results with continuum ChPT we have parametrized the chiral logarithms in terms of ij and 



£i, defined in eqs. (23) and (24). Introducing the scale independent coefficients Zj is useful in 
continuum ChPT, since the pion mass is constant in nature; consequently, the coefficients 
li are constants too. In lattice simulations we are free to choose the quark and pion mass, 
and it is this freedom that eventually enables us to compute Gasser-Leutwyler coefficients. 
For this it seems more useful to introduce the (scale dependent) coefficients 

fM 2 \ A 2 
h(fi) = l i + \n = ln^, (32) 





A 2 








fl 2 




^2 




= ln^| 




fi 2 



6(/i) = & + In I -jjj-J = ln^, (33) 

which are, up to irrelevant constants, the standard renormalized Gasser-Leutwyler coeffi- 
cients defined in [23]. 11 In terms of these coefficients the chiral logarithms are explicit and 
the results for the scattering lengths can be written in a more compact form. 

For simplicity let us first assume that we are interested in analyzing data at one fixed 
lattice spacing. In this case various terms which depend only on powers of the lattice spacing 
but not on the pion mass can be combined to single unknown parameters. This reduces the 
number of fit parameters. Explicitly we can write 12 



_ 7M 2 f Ml | ' Ml 40 \ 2c 2 a 2 [ 90 , M 2 

l ° ~ 327T/2 32ttV 2 I /i 2 21 ™ / + 32vr 2 / 2 I 7 % 2 



5 - 2C2aV «o 2 + ^|llln^M ; 04) 



32vr/ 2 V. 32vr 2 / 2 [ \i 



11 The constants are 7,/(327r 2 ). The are given in eq. (9.6) of Ref. [24] . 

12 The results for a\ are summarized in appendix |c[ 
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«21 



Ml ( Ml [1 Ml 4 lT=2 l 2c 2 a 2 f M 2 



3 In 



u 16tt/ 2 V 16vr 2 / 2 [2 /i 2 3 J 16tt 2 / 2 [ ^ 2 
2c 2 a 2 / 2c 2 a 2 fll, M 2 \\ 

The new parameters comprise the analytic terms through 0(a 2 ), hence these are con- 
stants for fixed a, except for the fact that they are scale dependent since they contain the 
parameters £»(//). In the limit a — > they assume the value = 1. As a short hand 
notation we introduced 

e° = ix + 2/ 2 - |/» + y , £ 2 = ix + 2/ 2 - jjl s + I , (36) 

for the combinations of LECs entering the scattering lengths. We emphasize again that we 
have expressed our results in terms of the decay constant in the chiral limit, and not in 



terms of /„-. If the latter had been used the LEC U [24] would also appear in eq. (36) 



Each of the expressions eqs. (34) and (35) contain five unknown parameters: the contin- 
uum parameters / and l l nn as well as kh,ki 2 ,c 2 , i.e. three more than the continuum result. 
This is already quite large, taking into account that one usually has data for a few pion 
masses only. In 0(a) improved theories we have Kn = 1 + 0(a 2 ), so in this case one may 
try to ignore the higher order corrections and set Kn — 1. 

We remark that the ratio Oq/M 2 has the functional form 

n l A lr, ]\/f 2 

M L 2 = W* +AW + A ™ Ml + A ™ Ml ^ Ml + A *° Ml + Ai °l^- (3?) 



The first two terms on the right hand side correspond to the tree level result in eq. (11). 



The constants A 00 — A 40 represent the five independent fit parameters, while A 40 is not 



independent. The first three terms in (37) have already been used in analyzing numerical 



lattice data jH], but the data could not be fitted well. It might be interesting to repeat 



the analysis with the full result in (37), even though the data was obtained for heavy pion 
masses between 500 MeV and 1.1 GeV, and ChPT is not expected to be applicable. 

If one wants to simultaneously analyze data for various lattice spacings one has to keep 
in mind that the coefficients kjj are no longer constants but functions of the lattice spacing, 

k v = 1 + K^a + K,[fa 2 . (38) 

This increases the number of free parameters in each formula from five to seven for unim- 
proved theories. In 0(a) improved theories we have k$ = and the number of free param- 
eters is increased by only one. Note that for unimproved theories the 0(a) corrections are 
not independent, 

7Kqi = — 5/?2l\ ^02 = K 22 ■ (39) 

Hence, if data for both scattering lengths are available the parameters in a simultaneous fit 
are related. 
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The fit formulae for the GSM regime are obtained from eqs. (34), (35) by dropping the 
appropriate higher order terms, leading to 



32tt/ 2 

Ml 
16tt/ 2 



K 2l 



Ml 



327T 2 / 2 

Ml 

16tT 2 / 2 



M 2 
51n^ 


40 
~ 21 


2 /i 2 


-1' 



5 • 2c 2 a 2 
' 32vr/ 2 

2c 2 a 2 
16tt/ 2 ' 



(40) 
(41) 



Here «/y = 1 + (tjjV The free fit parameters are /, c 2 and Ky , and the latter vanish in 0(a) 
improved theories. These formulae lead to the general form 



Ml 



^ + A w + A 20 Ml + A m Ml In M%, 



(42) 



for the ratio af,/M 2 , in contrast to (37). 



IV. CONCLUDING REMARKS 

Present day lattice simulations are still done with quark masses much heavier than in 
nature. Therefore, a chiral extrapolation to the physical point is still necessary, and that is 
where the predictions of ChPT enter the analysis of numerical lattice data. However, the 
results of continuum ChPT get modified at non-zero lattice spacing. For pion scattering 
with Wilson fermions we essentially find two modifications. First, the I = 0,2 scattering 
lengths do not vanish in the chiral limit, but rather assume a nonzero value of order a 2 . 
Second, additional chiral logarithms proportional to a 2 appear in the one-loop results for 
these quantities. Ignoring these modifications and using the results of continuum ChPT 
is potentially dangerous, depending on the size of these extra contributions. One either 
introduces a systematic error in the chiral extrapolation or the data cannot be fitted at all 
with the continuum results. 

Related problems may arise in the determination of the Gasser-Leutwyler coefficients 
associated with pion scattering. Experience shared by many lattice groups is that the lattice 
data does not show the characteristic curvature due to the continuum chiral logarithms. A 
possible explanation is the presence of additional chiral logarithms proportional to the lattice 
spacing. These can conspire with the continuum chiral logarithms such that the overall 
curvature of the data is diminished. A correct and precise determination of the Gasser- 
Leutwyler coefficients using the continuum results is very unlikely in that case. Using the 
expression derived here should help in that respect. 

The formulae we presented in this paper involve the parameter c 2 , which also determines 
the phase diagram of the theory. In particular the 1 = 2 scattering length may provide a 
handle to obtain an estimate of c 2 . At least the sign of c 2 should be easily accessible, and 
this is what matters for the phase diagram. 

We conclude with a remark on finite volume corrections. The numerical calculation 
of phase shifts and scattering lengths is usually done by employing the so-called Liischer 
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formula [4*2l |4"3] in order to circumvent the Maiani- Testa no-go theorem jH] . This formula 
relates the two-pion energy eigenvalues in finite volume to the infinite volume scattering 
length of the two-pion scattering process. 

In addition to the power law finite volume dependence of the the two-pion energy eigen- 
values, which one exploits to extract the infinite volume scattering lengths, there are expo- 
nentially suppressed finite volume corrections. These have been studied in Ref. [45] and it is 
straightforward to include these in our results. However, these corrections are expected to 
be very small on typical lattice sizes, much smaller than the corrections due to the non-zero 
lattice spacing. For example, for a pion mass of approximately 300 MeV and a finite volume 
with L ~ 2.5 fm the finite volume correction to the 1 = 2 scattering length is about 1 to 
2% [45J. 13 On the other hand, taking \2c 2 a 2 \ « (185 MeV) 2 at a « 0.086 fm found by the 
ETM collaboration at a charged pion mass of about 300 MeV [32J |33] , we obtain the rough 
estimate of about 35% for the 0(a 2 ) corrections in this channel. Although the numerical 
value for 2c 2 a 2 has a large error bar and can easily be a factor of 2 or 4 smaller, these 
numbers indicate that the finite volume corrections are most likely much smaller than the 
lattice spacing corrections. 
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APPENDIX A: HIGHER ORDER TERMS IN THE CHIRAL LAGRANGIAN 



The terms required for one loop calculations in the LCE regime are shown in (14). Among 
them are the NLO terms of continuum ChPT, which are given in Ref. [21]. The terms of 
0(p 2 a, ma, a 2 ) can be found in Refs. [15] 122]. So far unknown are the contributions of order 
p 2 a 2 , ma 2 , a 3 , a 4 . However, these missing terms are easily constructed with the spurion fields 
introduced in [2"2"] . 

There are essentially two independent spurion fields, M and A. These stem from the two 
sources of explicit chiral symmetry breaking, the quark mass and the lattice spacing. Under 
chiral symmetry transformations these fields transform as 

M — > LMR\ A — > LAR ] . (Al) 

These fields together with the field S and its derivatives are used to write down the most 
general chiral Lagrangian that is compatible with chiral symmetry, parity and charge con- 
jugation. Once the terms in the chiral Lagrangian have been found the spurion fields are 



13 For this estimate we made the crude approximation k cot 5 ~ 1 /a} 
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assigned to their physical values, 



M — > diag(m u ,m d ) = ml, A — > al , (A2) 

where / denotes the two dimensional unit matrix (recall that we ignore isospin violation 
and assume m u = ma = m). The number of terms in the chiral Lagrangian can be reduced 
further by making use of the Cayley-Hamilton theorem. 

Proceeding along these lines we find the following independent terms of 0(p 2 a 2 ,ma 2 ): 

C [p 2 a2] = ai a 2 (d^d^} + a 2 a 2 (^£^ t £t)(£ + £t)2 

+ a 3 a 2 (^(S + St))(^(E + £*)) , (A3) 
£ [ma2] = b 1 ma 2 {E + ^) + 6 2 ma 2 (E + E t ) 3 . (A4) 

The coefficients (Xj , bj 3X6 undetermined LECs. 

Since the spurion fields M and A transform identically under all symmetries we trivially 
obtain the 0(a 3 ) terms by replacing M with A in the 0(ma 2 ) terms. This leads to 

£ [a 3] = dxa^E + Et) + d 2 a 3 (E + E*) 3 . (A5) 

Finally, for the 0(a 4 ) terms we obtain 

C W = ei a 4 (E + Et) 2 + e 2 a 4 (E + E f ) 4 . (A6) 

Here we have dropped the constant term eoa 4 , which is also admitted by the symmetries 
but does not contribute to the pion mass and scattering amplitude. 

Note that we have obtained the 0(a 3 ,a 4 ) terms although the Symanzik effective action 
has not been constructed to these orders. However, any spurion field at these orders has to 
transform as an appropriate tensor product of the spurion field A and therefore gives rise to 
the terms in (A5), (A6) only. 



APPENDIX B: THE SCATTERING AMPLITUDE TO ONE LOOP 



The form of the one-loop result for the scattering amplitude is defined in eq. (20). The 



functions B(s, t, u), C(s, t, u), introduced in eqs. (21 ), (22), contain the one-loop corrections. 
The parts -B C ont ; Ccont ar e the contributions from continuum ChPT and read 



B 



cont 



96tt 2 / 4 



a 



cont 



3(s 2 -M 4 )F( S ) 

+ {t(t-u)- 2M 2 (t - 2u) - 2M 4 }F(t) 

+ {u(u-t)- 2M 2 (u - 2t) - 2M 4 )}F (u) } , (Bl) 
967r2/4 1 2 + 2/3) (s - 2 Ml) 2 + (/ a + 7/6) (s 2 + (t - uf) + M 4 ^ 
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The constants U are defined in eq. (23), Mq is the tree-level pion mass of eq. (J5J), and the 
function F(x) is given by 



Fix) 



-a In 



1 + a 



a 



a 



1 



4M 2 

x 



(B2) 



These results agree with the ones in [39]. The only difference is the use of the function J(x) 
in [32], which is related to our F(x) by 



16tt 2 J(x) 



Fix) + 2. 



This leads to some differences between our functions B cont , C cont and the ones in 
result for the scattering amplitude A(s,t,u), however, is the same. 
For the 0(a 2 ) contributions we find 



(B3) 
The 



B a2 (s,t,u) 



1 



96vr 2 / 4 



C a 2(s,t,u) 



967T 2 / 4 



(42c 2 a 2 + 18M 2 - 2As)F (s) 
+ (12c 2 a 2 - 12M 2 + 6t)F(t) 
+ (12c 2 a 2 - 12M 2 + 6u)F(u)} , 
- 30s(& + 1) + 3M 2 (11£ 5 + 6) 



+ 66c 2 a 2 (£ 6 + 1) + k 



(B4) 



(B5) 



The constants = ^(Mq) are defined in eq. (24) and fcs is the LEC associated with the 
0(a 3 ) terms (it is a combination of the Cj in (A5)). As claimed in section III C these 
functions assume a finite value for a — > 0. Hence the scattering amplitude reduces to the 
continuum result in this limit (recall the additional factor 2c 2 a 2 in front of B a 2,C a 2 in eqs. 

©, dH). 



APPENDIX C: THE 1 = 1 SCATTERING LENGTH 

With the result for the scattering amplitude of the previous section we can compute other 
quantities like the slope parameters bj or the phase shifts 5f . Here we only quote the one- 
loop result for the scattering length a\ in the isospin one channel, because the corrections 
due to a non-zero lattice spacing are fewer than in the other two isospin channels. 

Performing the partial wave expansion for the 1 = 1 case we straightforwardly obtain 

•i - & - H> - «*> + S} - - f }) 

^^(^ («) 



24?r/ 2 V16tt 2 / 2 [12 
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for the scattering length. We recover the continuum result in Ref. [23j for a vanishing lattice 
spacing. 

Note that the number of lattice spacing corrections is reduced compared to the results 



(25), (26) for the other two isospin channels. There are no corrections proportional to a 
and a 3 in the result for a\. Also the a 4 In M5 2 term is missing. The absence of these terms is 
a direct consequence of taking the difference A(t, s, u) — A(u, t, s) for the 1 = 1 scattering 
amplitude. However, analytic corrections of 0(a 3 ) enter when we replace Mq by M 2 and 
express the scattering length as a function of the pion mass: 

M - <- M - ^UmD-UmD-IumD^ 



24tt/ 2 V 12tt 2 / 2 [ 1V zv 8 ° v WJ 48 

, W a 2c 2 a 2 [ 5-. ~. 35 



f 2 16tt 2 / 2 [ 2^ oy *' *' 6 



24tt/ 2 V P P lQn 2 P [12 



For completeness we also give the analogue of the fit formulae in section HI E Replacing 



the U and £3 using eqs. (23), (24) we can write 



_ mi f_ mi is it; „, I=1 i 2^ r 15 . Ml 



2c 2 a 2 



«42, (C3) 



24tt/ 2 

where we introduced the combination 

The coefficients n±j comprise the analytic terms through 0(a 2 ). They are of the form 

n n = 1 + K^a + K^a 2 , K12 = n^a + nf^a 2 ■ (C5) 

These coefficients are constants for a fixed lattice spacing. The coefficient is related to 
the corresponding coefficients for the other two isospin channels: 

7 '^oi^ = ^^21 = 5ac^^ . 

In the 0(a) improved theory all these coefficients vanish. 

There is no relation for k^ 2 since it does not involve the low-energy constant ^5. Note 
that K12 vanishes for a —>■ 0, in contrast to the other two isospin channels. 

Finally, in order to obtain the result for the GSM regime one has to drop the a 2 M 2 InM 2 
term and sets Kn = 1 and K12 = 0. The resulting expression is simply the one-loop result of 
continuum ChPT. 
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